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Abstract 

We investigate the radiative decay of the axion into two photons 
in an external electromagnetic field to one loop order. Our approach 
is based on the world-line formalism, which is very suitable to take 
into account the external field to all orders. Afterwards we discuss 
how the calculation could be generalized to finite temperature. 



1 Introduction 



In order to solve the strong CP problem Peccei and Quinn [1] proposed 
that the full lagrangian of the extended standard model is invariant under a 
global chiral U{1) symmetry, which however is spontaneously broken, leading 
to a pseudoscalar Goldstone boson [2], the axion. This particle couples to 
fermions, depending on the model either to the known fermions of the stan- 
dard model (DFSZ model [3]) or to some exotic heavy new fermions (KSVZ 
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model [4]). In the latter case the axions couple to ordinary matter and radi- 
ation mainly by means of the triangle graph (fig. |l|), leading to an efi^ective 
two-photon (or two-gluon) coupling. 




Figure 1: Triangle graph 

Although the original Peccei-Quinn model, which assumed a symmetry 
breaking at the electroweak scale, is 'ruled out' by experiment [5], models 
with a much larger breaking scale, leading to light and weakly coupled, so 
called invisible, axions are still of great interest. Their astrophysical impact, 
especially their infiuence on stellar evolution, is discussed e.g. in [6-8]. 

In the present paper we study the axion-two-photon coupling in the pres- 
ence of a constant homogeneous electromagnetic background field. Thus we 
assume a pseudoscalar coupling of the axion field 0' to fermions, \'(j)'ip'y5ip, 
although the Goldstone character of the axion would in general require a 
derivative coupling. It is argued, however, in [6] and references therein, that 
the pseudoscalar coupling leads to the correct results as long as only one 
axion is attached to a fermion line. This is obviously the case for the triangle 
graph. 

Our calculation is based on the world-line formalism [9-17]. This is a first 
quantized formalism, originally used for one loop calculations in quantum 
field theory. The generalisation to higher loop orders is however possible 
as elaborated in the scalar theory and particularly in QED (see ref. [14- 
16]). Although the idea of dealing with particles on a world-line is quite old 
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[9], it has experienced a renaissance through the work of Bern and Kosower 
[10], who derived a set of new Feynman rules for one loop gluon scattering 
amplitudes from string theory, which thus inherited its superior organisation 
of perturbative calculations. Strasslers observation that the new set of rules 
can alternatively be derived through pathintegral methods has been very 
useful [11]. In this approach one writes the one loop effective action as a 
(super)particle path integral and evaluates this path integral using world- 
line Green functions appropriate to a one dimensional field theory on the 
circle. 

The world-line method is particularly useful in situations with a constant 
background field [16] and has been successfully applied to the calculation of 
photon splitting in a strong magnetic field [18]. The description of scalar 
and pseudoscalar couphngs on the world-hne is given in [19-21]. In [20] the 
world-line lagrangian was derived by analogy to the expression in the second 
order description of the fermionic one loop effective action in a background 
of vector, axial vector, scalar and pseudoscalar fields. Thus it was mathe- 
matically necessary to introduce two further bosonic and fermionic degrees 
of freedom. This procedure was clarified in [21], where the world-hne action 
for a fermion coupled to vector, scalar and pseudoscalar fields is derived from 
an eight dimensional Dirac operator in six spacetime dimensions coupled to 
a gauge field whose first four components are those of the vector gauge field 
while the fifth and sixth components are the scalar respectively pseudoscalar 
field. The Dirac fermion however only depends on the first four coordinates. 
This interpretation is possible if one doubles the fermion degrees (without 
altering the value of the effective action) and introduces the six Gamma 
matrices 

which obey the usual anticommutation relations 

{^A^'^b} = '^^ABhxS (2) 

The above mentioned six spacetime dimensional Diracoperator then has the 
form: 

E - r^(p^ - eA^) - r,X<p - TeAV (3) 
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Due to the presence of the pseudoscalar field, the effective action has both a 
real and an imaginary part, depending on whether an even or an odd number 
of pseudoscalar particles are involved. The real part of the effective action 
can be calculated as 

TsR = --InDetp] = - T — f Vx f D^e"^ (4) 
2 Jo S JP JA 

where P and A denote periodic and antiperiodic boundary conditions resp. 

and S is the world-line action, which for the case at hand (i.e. in the formulas 

of [20] we discard the axial vector contribution, introduce the fermion mass by 

setting the scalar field to — m/A and integrate out the auxiliary variables 

and xq via their 'equation of motion') is 

S = rrfr{^ + ^7/;-^ + ^^5V^5 + ^V^64 + m2 + A'20'2-2iAV6V^^<9^0' 
Jo 4 z Z Z 

+ie{x^A^ + 2i)^i,,d,A^)} (5) 

As the triangle graph involves only one axion field it contributes to the imag- 
inary part of the effective action. According to [20] its derivative with respect 
to the pseudoscalar field is given byQ 

j-jtT^ = - — Vx V^pn^^e-^ (6) 
dcp' Jo s Jp Jp 

where S is again given by (|^), the boundary conditions for the fermions have 
changed to periodic ones (as a consequence of the appearance of 75 in the 
usual Feynman calculation) and the insertion operator is (again for the case 
at hand) 

Q^, = -iX' dr Qv^x^Ve - imip^ipe^ (7) 

When solving the fermionic pathintegral by means of Wick-contractions with 
the world-line Green functions, we have to take into account the periodic 
boundary conditions of the fermions. There are now fermionic zero modes 
and the world-line supersymmetry is no longer broken by the boundary con- 
ditions. Thus we have 

^The normalization stems from a calculation of the anomaly in four dimensions (see 
e.g. [32], p. 552) and agrees with the one of [21]. 
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G^;^'"'''' = Gb , Gb{t„T2) = \t,-t,\-^^^^—^ (8) 

s 

where Gp/s means the fermionic respectively bosonic Green function on the 
world-hne. 

2 The case {F^^ ^0, T, = 0} 

The amphtude for the couphng of the axion to the electromagnetic field is 
given in the world-line formalism by: 

i-iX') dr' Q^;.y^^6 - imV^sV'e) e^^y^^'^e^^^^'^-^"' (9) 
exp \~j dr 



ip' \ • 1 • 1 • 



where the integration over the zeromodes has been seperated (a; = x^^^ + 1/, 

and l^dTy^ir) = 0). One has to use ■?/'aj,5,6 = "^^.5-6 + "^jt^s.e the integrand. 
As the integration over the fermionic zeromodes is only nonvanishing if the 
integrand contains a factor ■il)^^^il)':^^il)^^ip^^^il)f''il)f\ where the greek indices 
run from 1 to 4, producing an e^yapi one realizes, that the first term of 
the inserted operator does in fact not contribute. One also notices, that 

dr '0^°5,6'0/i,5,6 = so that in the kinetic part of the fermionic world- line 

action the substitution ip^^^^ '0/i,5,6 is allowed. 

As described in [16] we can seperate the electromagnetic field into the 
constant, homogeneous background field F and the outgoing photon fields F 
(which are described like the axion by plane waves): 

F^, = F^, + F^, (10) 

Using the Fock-Schwinger-gauge for the constant field yields: 

A,{t) = 4 + 4(t) = ^y^F,^ + e^,e'^-<^^ + e^^e'^^^^^ (11) 



5 



We realize that the Fock-Schwinger-gauge is compatible with the Lorentz- 
gauge, so that we can assume ei ■ ki = 0, i = 1,2. The influence of the 
constant background field can partially be incorporated by a redefinition of 
the Green functions and the fluctuation determinants. In doing so we have to 
keep in mind the periodic boundary conditions of the fermions (see appendix 
A). 

If we expand the interaction part of the exponential in the integrand 
of (P) and keep only terms, which are linear in ei and €2 and which do 
contribute after the integration of the fermionic zeromodes, we are left with 
the computation of the following expression: 



Jo 

1 



exp 



dr 



y2 1 _ 1 . ^ . _ . 



(12) 



where we have used / dr if)'^^'^ F^y^^ir) = and / Vtp^/e exp{— Jq dr IV's/e'^s/e) 
= 1. The rs are integrated over the circle with circumference s and one r 
can be set to zero for convenience because of translational invariance on the 
world-line (in our calculations this will be t^). We now perform the Wick- 
contractions with the bosonic Green function^ in a constant background field, 
i.e. we apply the following rules: 



< 2/M(n)y,.(r2) > 

< yt,{ri)y,y{T2) > 



-G^iuin, T2) 

-Q'^,An,T2) 
-^^M-(ri,r2)A;,e^^^(^^) 



(13) 
(14) 
(15) 



In the following Q always means the bosonic Green function Qb ( A. 2 ) 
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^^^kMr^)^^kMr2)^^kMr^)^ ^ exp{ ^ A;„^A;„,^;^, (r„, T^) } (16) 

m<n 

<M^i)Mr2)> = ^Qi..{ri,T2) (17) 

where the dot and prime mean differentiation with respect to the first respec- 
tively second argument. We get rid of the term G',j,„{ti, T2) — \ {pTiG'^xvi'^i-i '^2) 
-\-dj^Q hi,{ti, T2)) by partial integration, using the symmetry properties: 



Qh.u{Ti,T2) = -Qvn{T2,T-i) (18) 

QA-ri.r^) = -6^;.(ti,T2) (19) 

and noticing that the boundary terms vanish because of the periodicity of 
^/ii^(0, T2) = Qij,u{s, T2). The unbroken supersymmetry on the world- line leads 
to extensive cancellations in the calculation. We finally get as our result 
(setting k^ — q for convenience): 

M = Mo + Mi + M2 (20) 

with 



= AmX'e^ — dTidr2dT3e-'^'{27ry5{ki + k2 + k3){A7rs)-^ 

Jo s Jo 

GXp{— ^ ki ■ Qiji, Tj) ■ %}e;uii/i;U2t'2^1/ii^li'i^2^2^2fc'2 (21) 

fCO fS 

Ml = 2imX'e^ ds dndT2dT3e-'^ '{2n)^d{ki + k2 + k3){4ns)-^ 
Jo Jo 

2{ki • F • ei) {(62 • g{r2, n) ■ ki) + (e2 • ^(t2, T3) • A;3)} 
-{ki- F ■ k2){e-i ■ Q{ri,T2) ■ 62) - (ei • F • 62) (/ci • Q{tuT2) ■ /C2) 



+ 2{ki-F-e2)iei-g{n,T2)-k2] 
+ (ei, ki) ^ (£2, /C2) 



(22) 
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"1 n n ^ 

M2 = — mX'e^ dss dTidT2dT3e'"'' '{27c)^ 6 + k2 + h){A7Ts)-^ 
2 JO JO 

^(ei ■ G{ri,T2) ■ £2) [{ki ■ G{ti,T3) ■ k^) - (/ca ■ G{t2,T3) ■ ^3)} 
+ 2{e, ■ g{n,rs) ■ k:,){e2 ■ g{r2,n) ■ k,) 

+ (ei • g{n,T3) ■ k^){e2 ■ g{r2,n) ■ k^)] (23) 
+ (ei, ki) (62, /C2) 

where we have introduced the notation Ffj^^ = ^e^^a/sF'a/s- 

Before going on to the finite temperature case we want to specify our 
formula to the case described in [22], i.e. we neglect the axion mass and 
assume the crossed field case {E ± B,E = B or F^^F^^, = F^^F^y = 0). In 
the limit of a vanishing axion mass we additionally have q || ki || k2 because 
of energy- momentum conservation, that means ki = —Aq, ^2 = —(1 — -A)q. 

We see immediately that under these circumstances only Aii contributes 
to the amplitude. Furthermore we easily verify, that F^aFapFpi, = 0, so that 
we can use the expansions (|A.5|) and (|A.6|) in A^i. In doing so we arrive at 
the following formula (see appendix B for some details): 



A 



M = AimX'e^ ds / rfridr2C^r3 e-'"''(27r)^(5(fci + ^2 + fc3)(47rs)-2 

JO JO 1 — 

^^P{^H^i ■ ^(^i'^i) • kj}{ki ■ F ■ ei)(e2 • [^(r2,ri) - ^(r2,r3)] • k2) 

+ {ei,k,)^{e2,k2) (24) 

= -mX'e\-\2n)^5{ki + k2 + q) ds dai das e""^ ^ 

4 JO JO JO 

Y^^{fifiuFf,^){f2fiuF^u){(^i - 2cri(T2 - (rl) + (ei, ki) ^ (e2, ^2 
exp{je\q ■ F ■ F ■ q) [-{a, - alfA' - (^2 - alf{l - Af 

2(Ti(T2(l - 3(Ti - 3(T2 + 2(T? + 2(T^ + 3(Ti(T2)A(l - A)]} (25) 

with 
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fiaf3 — kia^ifj Giaki/3 (26) 

Equation ( pS]) can be shown to be equivalent to the result obtained in [22], 
up to an unphysical overall minus sign. To see the equivalence of the two 
formulas one would have to transform the variable s according to 

s s{e2(g-F-F-g) [2(TiCT2(l-3ai -3a2 + 2a^ + 2a^ + 3(Ti(T2)A(l-A) 
- {a, - alfA' - {a, - a^l - Af]}'/' (27) 

perform a Wick-rotation in s and change from Euclidian to Minkowski space. 

3 Remarks about the case {F^^jj 7^ 0, T, /i 7^ 0} 

The finite temperature [25,26,33] and density case even without background 
field raises some new questions and we first treat this case. Finite temper- 
ature is usually taken into account on the world-line in the context of the 
imaginary time formalism. Thus the time component becomes a circle with 
circumference (3 = ^. Under this topological constraint the pathintegral in 
(P), which is over closed paths, can be written as a sum over pathintegrals 
without the topological constraint but each one over paths where the end- 
point differs from the starting point by a multiple of (3 in the time component 
[23,24]. This takes into account the different windings of the closed paths 
on the time-circle. In order to get the statistics right the sum has to be 
alternating^. That means, if P^y is the pathintegral over paths from y to x 
under the topological constraint = + nj3 and is the one without 
the topological constraint, we have: 

oo 

Pxy = X! (~l)"-^(^,X4+n/3)y (28) 
n=— oo 

Thus we are led to calculate the following expression: 

■^Note that this is not in contradiction to the periodicity of '0(t) on the world-hne since 
ipir) is a world-hne and not a space-time fermion. 
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M = -4mA'e^e^^o,^ei^e2,./ciaA;27 / — X! / ^ Vx V^) (29) 



where we have already performed the integration over the fermionic zero 
modes. We now have to spht the path into a hnear and a periodic part and 
separate the center of mass 

x{t) + 71/3-64 + y{T) (30) 

s 

where 64 is a unit vector in the direction of the fourth (i.e. time) component. 
The remaining pathintegrals run over periodic paths with center of mass zero 
so that the usual Green functions can be used for the Wick contractions (this 
is true also for the fermionic pathintegrals). The result we get is: 



M = -4mA'e^e/,i.„^ei^e2,.A;i„/c27 (27r)'^5(A;i + A;2 + g) (5^,^+^,^+,^, 

/ dss^ daida2 V (-l)V('"'=i'^^+'"'=2<^2)"'^e-^e-'" *(47rs)-^ 
•^0 Jo 

exp{s(/ci • /C2GS12 + ki ■ qGbi + k2 ■ qGBi)} (31) 
where T3 has been set to zero again, 

Gsij = Wi - o-jl - (o-j - ajY , Gbi = cfi - af (32) 

and ujk are the external bosonic Matsubara frequencies. We now perform a 
Jakobi-transformation 



(-1)" exp{i7r[Tn^ + 2zn\} = (33) 

i=— 00 

= {—ir)~^^^e ~ ^2 exp •! ivr 



n=— 00 
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and get0: 



/■I °° poo _ _ 

/ daida2f3^^ / cis s^^^ exp{— s[m^ — /ci ■ A;2G'bi2 — fci ■ gC^i 
-'0 ^rr'^ Jo 

(34) 



- ^2 ■ qGB2 + (t^n - t^fciO"l - t^fcaCTa)^]} 



1 rl — <J2 

da2 / ' 
Jo 



dai exp{— s[m^ + kfai^l — ai] 



+ / d(Ji / da 2 eyi'l) {- s[m^ + kla lil - a i) + k^a 2(1 - cr 2) 
Jo Jo 



+ 2kik2cria2 + {uj^ - Uk^ai - Uk^il - 0-2))^]} 



(35) 



where the "— " in uj~ characterizes the Matsubara frequencies as fermionic 
ones. We can now introduce the chemical potential via a;^ ^ — ifi. We 
convert the sum in (^^ into contour integrals via: 



1 

27Ti 



ioo+/i— e 



2ni 



du —-, r / du ^, . 

■too+/i+e fi'-^) fioo+iM—e 



du ■ 



dujf{u) - i ^ du f{u) 



du 



g/3(/^-w) _|_ ]^ 



(36) 



For the case /i > the contour C' is depicted in (fig. |]) and we have 
assumed that fiu) is analytic in a neighbourhood of the line 'Sliuj) = /i and 
lim^^^-i-oo /(o- + ^-R) = for < a < /i (which is true for the case at hand). 

■^Because of the analytical properties of the integrand it is possible to change the order 
of the s-integration and the summation. This is also confirmed by formula (^2|) derived 
with the Matsubara formalism. 
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Figure 2: Contour C' 



In the following we consider only the case < m. Then our jiiS) is 
analytic inside the contour C and therefore the last term of ( PB]) does not 
contribute. 

We are thus led to calculate the following expression for the temperature 
dependent part: 



010 ^7 r 

2m J-ioo+p^+e e^^"^-/^) + 1 



1 |-ioo+^t-e 



iB) 



-3/2 



(37) 



with 



A = - ki- k2GBi2 - ki ■ qGei - k2 ■ qGB2 

B = ujk.ai + uJk^ai (38) 

We now have to deform the a;-contour to the real axis. In doing so we have 
to be careful because of the squareroot cut in the integrand. The details of 
the calculation can be found in appendix C and here we just give the result: 
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1 ^ 

= -^m\'e^e^^^^ei^e2ykiak2^{2TtY5{ki + k2 + q) 5^^^^ +0;^^ +c^, doida2 

dp ; 1 ; 



^-l3{-VA^/p2+l~fj,+iB) _^ ^ ^^-l3{-VA^/p^+l-ti+iB) _|_ ly 

(39) 

In appendix C we derive this result also with a different method. 

Decomposing (p^ + 1)^^ = —p^/ip^ + + 1 in the second and fourth 
term and performing a partial integration in the —p^jijy^ + 1)^^ part one 
arrives at a simpler expression^: 

M.^ = -^mX'e'^e^^a'yei^(^2ukiak2^i2Tc)'^6{ki + k2 + q)6^^^^+uJ^:^+^g daida2 



00 







dp PA 



-1/2 



^|3{^/A^yp^^+l-^l+^B) ^-f3(-VA^yp^2+l-^l+^B) 

+ 



^g/3(v^^pHl-M+iB) _^ ly ^^-P{-^/A^yp■2+l-fJ.+^B) _^ ^y 

(40) 



The general cx-integration is not possible here. A similar problem has 
been encountered during the investigation of a scalar quantum field theory 
at two loop level [31]. One could try to specialize to the case where the 
rest frame of the axion is identical with the frame of the heat bath and 
perform the analytical continuation in the external Matsubara frequencies 
{iukj (^i)o — i^j , j = 1, 2 , see e.g. [27] and for a general treatment 
of the relationship between imaginary-time and real-time finite temperature 
field theory also [28-30]). As pointed out in [26] the result should coincide 
with the corresponding calculation in the real-time formalism (see e.g. [25]), 
if we assume M << m, where M is the axion mass. The examination of the 
analytical continuation has not been finished yet. Some details can be found 
in appendix D. 



^This was suggested by H. Sato. 
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P, ai 



= r3 



p + ki, a2 



p-q, a-i 




Figure 3: Relationship between Feynman parameters and world-line times 



It is possible to show that our expression agrees with the one one gets with 
the usual Matsubara-Feynman-rules, if one introduces Schwinger parameters 
and integrates out the spatial momenta. In order to get (|^) one has to use 
the relationship between Feynman parameters and world-line times (see (fig. 



«! = Ti = sai , a2 = T2-T1 = s{a2-ai) , = S-T2 = s(l-a2) (41) 
In detail we have: 



M 



(2^ 
m — ^ 



tr 



m- (J + ^1 ) 



(c^„ + LOk^y + {p + kif + m? 



7m 



+ + m- 



r75- 



[uj- 



^qY + {p- q)"^ + 



-4:mef,ua-yeif,e2ukiak2'rP~^ V] / 77—7^ / daida2da3 

n=-oo-' {^T^r Jo 

exp{-ai{{u~f +]f + m^) - a2{{uj~ + ujk^f + {p + hf + m^) 
- "3((^;^ - ^^qf + {p- qf + m^)} 
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POO 

-4me^,,a^ei^e2z.fciafc27/?~^ / daida2da3 (47r(ai + 02 + "s))"^^^ 

n=-oo "^0 

exp{-m^(ai + ^2 + as) - (^^.« )^(ai + "2 + as) - ^fcia2 - ^^,^3 



, -/o o ^ ^ r2 I (^"3 - ^1^2) . 

+ (2u;g«3 - 2t^fcia2) - g as - A;ia2 H ■ ■ 1 

ai + a2 + as 

-Ameaua-yeiae2,ykiak2y(3~^ V / (is(47rs)"^/^s^ / da2 / dai 

exp{-s[m^ + {uj~ + uJk^{l - ai) + Uk^il - 0-2))^ + klai{l - ai) 
+ kla2il - a2) + 2kik2ai{l - a2)]} (42) 

where we have used tr(757^7,^ 70,7/3) = 4:e^i,ai3 with the Euchdian conventions 
ei2S4 = 1, 75 = — 7i727s74, 74 = "^70- In order to get the first term in the 
square brackets of (|35|) one has to transform ct2 — ^ 1 — o"2 and n — ^ — n in (|4^) . 
The second term in the square brackets of (^) is derived by interchange of 
Ti and T2. 

With the same technique it is possible to treat the case with a constant 
background field. The formula we derive along the lines outlined above is: 







X'tti -> -> /'°° 2 

+ujk2+^q I '^'^ d(7ida2dase * 

00 

exp{A;i ■ g{Ti, T2) ■ k2 + ki- Q^n, n) ■ q + k2 ■ Q{t2, n) ■ q} 

n=—oo 

exp{-s(a;~ - cUfciCXi - ujk2(^2 - t^gO-g)^} {A<o + Mi + M2} (43) 

with 

Mo = 2eh^'\h,uf2,u) (44) 
Ml = 4^eV/2[^(/2^,F/..)((ei-^(ri,r2)-A;2) + (ei-^(ri,r3)-g)^ 

-]^{ki-F-k2){ti-G{n,r2)-e2) 
+ {ki-F-t2){ti-G{n,T2)-k2) 

+ (^n - ^fciO"! - ^^/fc2^2 - ^qCr:i){ei)i{f2^,uF^y) 
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~(ei ■ F ■ e2){ki ■ G{ti,T2) ■ k2 

€ S ^ iyF^yF^y 



+ {(^i,ki 



{^2,k2 



(45) 



ei ■ G{.T-i,T2) ■ t2){ki ■ Q{ti,t^) ■ q) 



+ (ei ■ ^(n,i"3) ■ q){^2 ■ ^(T-2,n) ■ ki) 

+ |-- + 2(t^^ - Uk.cri - uJk^a2 - tOqa^f^ (£1)4(62) 



2{uj^ - uJk^ai - uJk2(T2 - t^gCTs) <^^{ki)4{ei ■ Gin, T2) ■ 62) 
(ei)4 ((e2 • Q{t2, n) ■ ki) + {€2 ■ g{T2, T3) ■ g)) } 



1 



2 

+ (ei, ki 



(ei ■ Qin, Tg) ■ q){e2 ■ G{t2, Tg) • q) 
' (e2, k2) 



(46) 



where we have used (^6]) and 



^^iMua-yfa-r- Q{'^i,Tj) IS of couise Under- 
stood as Q{sai, saj). It is again possible to fix (73 = 0. (€4)4 and {ki)^ mean 
the fourth (i.e. timehke) component of the corresponding vector in brackets. 

We have seen that the world-hne formahsm provides a powerful tool in 
treating problems with a constant electromagnetic background field. The 
examination of temperature effects still suffers from the problem, that the 
occuring integrals are not solvable in general and that the result has to be 
continued analytically. These problems have also been discussed in [31], 
where it was shown however, that the world-line method offers an illustrative 
way of finding the divergent parts of the 0g-two-loop-amplitude at finite 
temperature. It is too early to draw a final conclusion about the usefulness 
of the world-line method in finite temperature field theory. It might also be 
possible to combine the real time formalism with the world-line method, so 
that no problems of analytical continuation occur. Here investigations are in 
progress. 
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A Green functions and determinants in a con- 
stant background field 

In a constant background the defining equation for the bosonic Green func- 
tion becomes 

\ (s,x^ - 2^eF,A^) G'nin, r,) = 8,. (^(n - r,) - (A.l) 

In the case at hand the fermions obey periodic boundary conditions. Thus 
we have [16]: 

eF \^sin(eFs) s) ^ ' ' 

where 

GBi2 = |ri-r2|-^^^^^ (A.4) 

s 

is the well known bosonic Green function on the circle without background 
field. 

The first few terms of the Taylor expansion in F for the Green functions 
are: 

Qbi2 = GBi2-\GBi2GBi2seF+'-Gl^2{eFf + 0{F^) (A.5) 

Gfi2 = Qbi2 = ^512 + 2i{GBi2 - ^)eF + '^GBi2GBi2s{eFy + 0{F^) 

(A.6) 

where in the expression for Qbi2 the constant coincidence limit has already 
been subtracted (which is possible because of momentum conservation). 
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The free pathintegrals in a constant background and with the fermions 
obeying periodic boundary conditions are: 



Vx exp j-;^ ^ dT{x^ -\- 2ieXfj,F^,^Xi,)^ = 

= detp^[-(9^]detp^[I-2ieF9-^] (A.7) 
(47rs)-2 

= det/[a]det/[I-2ieFa-^] (A.8) 
^ V ' 

1 

where the prime at the determinant reminds at the fact, that the zeromodes 
are excluded. 
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B Derivation of the formulas (|23) and ( [25l ) 



In the case ki \\ k2 \\ q we have, because of the antisymmetry of F and Q, 
that the terms proportional to ki-F-k2 and ki ■Q{ti, T2) ■ ^2 do not contribute 
(reahze, that the order of ri and T2 in the argument of Q can be changed 
because we are integrating both over the whole interval [0,s]; thus (|18D can be 
applied to see, that the term does not contribute). In order to derive ( P^ we 
take the first term in the square brackets of p^ ) as it stands and the last one 
with (ei, ki) (e2, k2). Using the fact that ki = —Aq and k2 = —(1 — A)q 
we immediately get: 



2{k, ■ F • ei) {(e2 ■ ^(ts, n) ■ k^) + (es ■ ^(ts, rg) ■ g)} 
+ 2(A;2-F-ei)(e2-^(ri,r2)-A;i) 

2A(1 - A)(g ■ F • ei) [e2 ■ Qin, T2) ■ q - e2 ■ <3{r2, n) ■ q 
2 



+ 



A 



(A:i ■ F ■ ei; 



£2 



- ^(r2,r3 ■ k 



(B.l) 



the first summand of ( |B.1| ) vanishes upon integration over the as explained 
above. 

To get from (^) to ( [25| ) we rescale the = suj and set = 0. Fur- 
thermore we realize that in the expansion of €2 ■ {p{sa2, sai) — Q{sa2, 0)^ ■ /c2 
only the second term of ( |A.6D contributes. The first one is proportional to 
the unit matrix and we have ■ fcj = because of the Lorentz gauge choice. 
The third one is odd under interchange of world-line times. Thus integration 
of CTi and a2 over the whole interval [0, 1] yields zero. As already stated in 



the text it is easy to verify that F^aFapFjs^ 



0. 



Moreover in the exponent of exp{| Si^j h ■ Qiji, tj) ■ kj} only the third 



term of (|A.5|) contributes in the limit of vanishing axion mass. In this case 



g = 0, so that the first term yields zero and the second term vanishes because 
of the antisymmetry of F. 

If we now split the range of integration into (ai G [0,1], (T2 G [0,cri]) 
and ((Ti G [0, l],o"2 G [(Xi, 1]) respectively and perform a change of variables 
(cTi — > 1 — 0"! in the first case and 0-2 ^ 1 — c"2 in the second case) we see 
that both ranges contribute equally and we end up with formula (p5|). 
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C Derivation of formula (|S9|) 



Our starting point is equation (pTl). We concentrate on the first of tlie two 
integrals, tlie second one is treated in an analogous way. The aim is to deform 
the cu-integration-contour to the real depicted in (fig. |^). 




Figure 4: Deformation of the a;-contour 

First we make a shift in uj by iB and rescale by \/~A ^ . This is possible 
because \/~A > m > for all values of o'i/2-Q Thus we get: 

^Explicitly we have: A = rri^ + 2ki ■ fc20'2(l — fi) + — crj) + fc|(cr2 — erf) + 
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^ [A - (to- iBf] 2 

3 3 

oo+i52 [l-a;21"2 .oo-i<52 [l-cj2l-2 



+ A Id^ i_/¥'g/3((l-5,exp{M)v^-M+^i?) + l ^^-^^ 

In the last integral we have introduced the variable z = 1 — uj = 5ie'^'^ . We 
have also taken the limit ^2 — ^ here. The range of integration comes from 
the fact, that we have taken the squareroot cut along the negative real axis. 
As we are interested in the limit 5i — > the only contribution from the last 
integral is 



The first two integrals in ( |C.1|) give the same contribution. If we have a pos- 
itive imaginary part for u the imaginary part of 1 — cj^ is negative. Therefore 
we have in the first integral (1— a;^)~^/^ = (e~*''(ci;^ — 1))"^''^ = —i{u^ — l)~^^'^ 
and a relative minus sign in the second one. If we introduce the variable 

= tu^ — 1 we get 



-2iA-^ 



POO 



doj 



■1]" 



= -2iA~^ 
-2tA-' 



^P{ujVA-ij,+iB) _|_ ]_ 
.oo -2/ 2 _^ ;^\-l/2 
/ dp ; 

^Pi-ZA-ii+iB) 



+ 1 



(C.3) 



dp 



(p2 ^ 1)^3/2 
g/3(v^^p2 + l-M+«S) 



02 > (Tl (cTi <-> 0-2, fcl ^ ^2)- 



_j_ 1 (^e/3(v^VP^+l-M+iB) _^ X)2 J 

a 2) for (Tl > (72 and a similar expression for 
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where we have omitted terms which vanish in the hmit 5i — 0. We see that 
the singular part cancels the one from ( |C.2| ) and we are left with the first 
two terms of ( P^D if we insert our result into (P7|). 

The second integral of (|37|) is treated similarly and leads to the remaining 
two integrals of (0). 



As already mentioned in the text it is possible to derive (|39D in a different 
way. Startingpoint is in which we make use of the following identity: 



/ dss^'^e-"' = / dss^ / dp, 
Jo Jo y/n Jo 

We get: 



TT JO 



dp(p2+a)-2 (C.4) 



1 poo 

daida2 / dp 
Jo 



duj -r, ^ 

2ni J-ioo+f^+e e^('^-^) + 1 



1 



2m 



(A + p'^-(uj-iByy^ 
duj 



(C.5) 



where A and B are again given by (0). We now apply the residue theorem 
in order to calculate the a;-integrals and use: 



Res. 



B±y/A+p2 



'{A + p^-{uj-iBy)-^' 

p±f3(u)~fj.) 



If we rescale p ^ we end up with 



+ 



(C.6) 
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D Some remarks about the analytical contin- 
uation 



If we try to carry out the program proposed after formula (^) and addition- 
ally use = and (fcj)o = ^ where M is the axion mass and set M = 
in order to get the leading order in — , the cr-integration becomes trivial and 
(H) reads (for /i = 0): 

Ztt ^ ^ 

POO 

/ dp 
Jo 

where Ai^ means the leading order finite temperature contribution. We now 
perform a partial integration of the first term in the square bracket: 



m 



X _|_ gm/3i/p2_,_i 



_2 i/^ -r -l; e 

+ pm 



(D.l) 



dp{p^ + ir/^^^±^= (D.2) 



.oo / 2p 1 ^ p e^'^V^ 
/ dparsmhlp) — — h mp—- -—r , 



If we use 2arsinh(p) = 2 ln(A/p2 _|_ x + p) = In ( J -^g arrive at: 



27r2 

dp 



P ^ /vg + T + A 1 ^j3_3^ 



(p2 + X)2 \^ + 1 - X + e'^^Vp^ 

X gm/SY^P^+l / y „ 

+ m/3^ — l + ln(A/p2 + X+p)- 



p2 + 1 _^ gm/3Y^p2+l-j2 \ ^ + 1 , 

This does not coincide with the results of [25,26]. The first term in the 
square bracket is exactly the known result but we get the additional second 
term of ( p.3| ). A closer look at the calculation in the Matsubara formalism 
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shows the following: If one continues the external Matsubara frequencies 
also in expressions exp{—iPuJk^^^}, which is indeed 1, one obtains exactly 
the additional term. Such continuations play a role because the angular 
integration produces a singularity to be cancelled by further factors 

~M^. This subtlety in the analytical continuation should be burried in 
our case in the (j-integrations, which correspond to the angular integration. 
Unfortunately they seem to be rather difficult to handle. 
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